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We discuss the dynamics of a quintessence model involving two coupled 
scalar fields. The model presents two types of solutions, namely solutions that 
correspond to eternal and transient acceleration of the universe. In both cases, 
we obtain values for the cosmological parameters that satisfy current obser- 
vational bounds as well as the nucleosynthesis constraint on the quintessence 
energy density. 
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1 Introduction 



Recent measurements of the Cosmic Microwave Background (CMB) anisotropics by var- 
ious experiments [|| suggest that the Universe is flat, with a total energy density that 
is close to the critical density. However, there is strong observational evidence that 
matter (baryonic plus dark) can only account for about one third of the total matter 
density. On the other hand. Type I Supernovae light-curves indicate that the Universe 
is accelerating at present 0. All these observations can be reconciled if one assumes 
that the the dynamics of the universe at present is dominated by a negative pressure 
dark component, the main candidates being a cosmological constant and dark energy 
or quintessence @, ||, H; the latter is characterized by a time- varying equation of state 
parameter, wq = p/p, approaching a present value wq < —0.6. 

A commom feature of all quintessence models presented sofar is that the asymptotic 
accelerating behaviour of the universe is driven by the dynamics of a single field. There 
are, however, several motivations for studying the case of coupled scalar fields. Firstly, 
if one envisages to extract the potential suitable for describing the universe dynamics 
from fundamental particle physics theories, it is most likely that an ensemble of coupled 
scalar fields (moduli, axions, chiral superfields, etc) will emerge, for instance, from the 
compactification process or from the localization of fields in the brane in multibrane 
models or from mechanisms for the cancellation of the cosmological constant (see e.g. 
0] and references therein). Furthermore, coupled scalar fields are invoked for various 
desirable features they exhibit, as in the so-called hybrid inflationary models ^ and 
in reheating models |T^, [TT[] . Finally, it has been recently pointed out that an eternally 
accelerating universe poses a challenge for string theory, at least in its present formulation, 
since asymptotic states are inconsistent with spacetimes that exhibit event horizons [|l^. 



13 1 . Moreover, it is argued that theories with a stable supersymmetric vacuum cannot 



relax into a zero-energy ground state if the accelerating dynamics is guided by a single 



scalar field |T^, |T3| . The main argument relies on the fact that , in a supersymmetric theory, 
one expects that the asymptotic behaviour of the superpotential is given by Wi^cj)) = 
Woe'""^^"^, which, in order to ensure the positivity of the 4-dimensional potential V"(0) = 
8|(90iyp — 12|iy^| implies that |a| > ^/6. However, this value is inconsistent with the 



requirement of an accelerated Universe at present \a\ = + ujqq) /2 < 1.5 0] as data 
suggest that ojqq < —0.6 |T^. The situation is different in the presence of fields that do 



not reach their minima asymptotically. Indeed, in this case, the asymptotic behaviour of 
the superpotential would be better described by the function 1^(0) = WQe~°"^^'^F{(f),il)), 
where -F(0, ip) is a polynomial in the fields (f) and if) and the positivity condition then 
becomes: — 6 + 4[((90F)^ + {d^F)"^]/ F"^ > 0. One can then easily see that, by a suitable 
choice of the polynomial F{(j),ilj), the positivity condition can be reconciled with the 
requirement of successful quintessence. 

Recently, a two-field model has been proposed [jl5[, exhibiting a class of solutions 
where fields do not necessarily settle in their minima at present, thus evading some of the 



1 



conclusions of Refs. [|l^, jTSf , concerning the stability of the supersymmetric vacua. The 
potential is given by 



l^((^,^) = e-^^P(0,V') , (1) 

where 

P(0, ^) = a + {ct>- 0o)' + h{il,- i,o? + c - V^o)' + di,{ct>- 0o)' , (2) 

in units where M = (SvrG)^^/^ = h = c = 1. Such potentials arise in the low-energy 
limit of fundamental particle physics theories such as string/M-theory and brane-world 
constructions. The overall negative exponential term in signals that this could be moduli 
type field which has acquired an interacting potential with the ip field. 



2 Cosmological Solutions 

We consider a spatially-flat Friedmann-Robertson- Walker (FRW) universe containing a 
perfect fluid with barotropic equation of state P7 = (7 — 1)P7, where 7 is a constant, 
< 7 < 2 (for radiation 7 = 4/3 and for dust 7 = 1) and two coupled scalar fields with 
potential given by Eq. (0). The evolution equations for a spatially-fiat FRW model with 
Hubble parameter H = a/a are 

= -3H{p^+p^) , (4) 
= -3H^-d^V, (5) 
^ = -3H7p-d^V, (6) 

where 9<^(^)V" = subject to the Friedmann constraint 



H' = l[p, + l<P' + l4^' + V) , (7) 



The total energy density of the homogeneous scalar fields is given by pq = 0^/2 + iIj'^/2 + 
We define the dimensionless variables 



^ = ; y 



V _ i) 



The evolution equations can then be written as an autonomous dynamical system: 
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3 3 

x' = -3x - d-d^\n{V)y^ + -xf{x,y,z) , (9) 

y' = ^y{xd^HV) + zd^HV)) + ^yfix,y,z) , (10) 

z' = -3z-^d^\n{V)y' + ^zf{x,y,z) , (11) 
0' = V6x, (12) 

i^' = Ve z , (13) 

where f{x, y, z) = 7(1 — x^ — y"^ — z"^) + 2(a;^ + z^) and a prime denotes a derivative with 
respect to the logarithm of the scale factor = log(a/ao), where is the scale factor at 
present. The constraint equation becomes 

+x^ + y^ + z^ = l. (14) 



2>H^ 
and, therefore 

^Q^{§-,=x^ + y' + z' . (15) 

Integrating the above field equations, one finds that there are essentially two realistic 
types of behaviour, illustrated in Figure (notice that N=-30 corresponds to the Planck 
epoch; nucleosynthesis occurs around N=-10, radiation to matter transition around N=-4 
and N=0 today). Model I (A = 9.5, a = 0.02, 0o = 29,V^o = 15, 6 = 0.002, c = 6x 10-^ c/ = 
4.5) corresponds to the case where vacuum domination, which occurs when VLqq > 1/2, is 
permanent and Model II (A = 9.5, a = 0.1, 0o = 29, ^0 = 20, 6 = 0.001, c = 8 x 10-^ d = 
2.8) to the case where vacuum domination is transient. Permanent and transient vacuum 



domination have also been found in the model of Refs. |]T6[ and |jT^, respectively. Of 
course, there remains the (non-realistic) case where accelerated expansion never occurs. 

The effective equation of state for quintessence is given by pq = ujqPq, where wq = 
IQ - 1, and 



PQ v + (^72 + ipy2 

x^ + y^ + z^ ^ ' 

A necessary and sufficient condition for the universe to accelerate is that the deceler- 
ation parameter, q, given by 
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Figure 1: Evolution of Qq (full curve), Qr (dashed curve) and flm (dotted-dashed curve) 
for Model I (upper panel), corresponding to permanent vacuum domination and Model 
II (lower panel), corresponding to temporary vacuum domination. 



q = -^ = -{l + 3wQnQ + nr) , (17) 

where Qr is the fractional radiation energy density, is negative. 

In Fig. 1^, we show the evolution of wq and q Model I (uppper panel) and Model II 
(lower panel). 

In both models, the equation of state has reached wq ~ —1 for the present time, as 
favored by the available data |jl8| (and making it hard to distinguish from a cosmological 
constant) but, whereas in Model I uq remains negative, in Model II it is in the process 
of increasing today towards positive values, then oscillates slightly until it reaches its 
asymptotic value. Similarly, in both models, the deceleration parameter is negative today 
but, whereas for Model I q remains negative, in Model II it oscillates and becomes positive 
before it reaches its asymptotic value. 

Scaling behaviour is apparent in Figure ^, where we plot quintessence energy density 
Pq as a function of the scale factor for Model I, for different initial conditions, namely, 
Pq Pr, Pq > Pr and pq pr- After showing some initial transient, each solution scales 
with the dominant matter component before pq begins to dominate. Model II presents 
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Figure 2: Evolution of the equation of state parameter wq (full curve) and the de- 
celeration parameter q (dashed curve) for Model I (upper panel) and Model II (lower 
panel). 



similar behaviour but, as the scalar fields never stop rolling, matter-dominated scaling 
evolution is soon resumed. Our models seem to be more sensitive to changes in the initial 
conditions than models with just one scalar field (this is to be expected since there is 
more freedom e.g. in the way kinetic energy is shared between the two fields) but no fine 
tuning of the initial conditions is needed. 

Permanent vacuum domination takes place when at least the (p field ends up settling 
at the minimum of the potential, thus corresponding to a cosmo logical constant; in Model 
I, both fields settle at the minimum of the potential, see Figure ^. Transient vacuum 
domination occurs either when the potential has no local minimum or arrives at the 
local minimum with enough kinetic energy to roll over the barrier and resume descending 
the potential. Notice that the evolution of ip is slight compared with 0, especially for 
Model II. 

Both models satisfy present bounds on relevant cosmological observables. The tightest 



bound comes from nucleosynthesis, ^q{N ~ —10) < 0.044, requiring A > 9 [19| from 
the most recent CMB data, Qq < 0.39 at last scattering, is less stringent than the 
nucleosynthesis bound. Other bounds we take into account are: flm = 0.3 ± 0.1, wq < 
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Figure 3: Evolution of the background (matter -|-radiation) energy density (dashed) and 
quintessence energy density, pg, for Model I, for different initial conditions: pq -C pr 
(dotted), Pq ~ pr (solid) and pq > pr (dotted-dashed). 

-0.6, VLq = 0.65 ± 0.05 and h = 0.65 ± 0.05 H today. Indeed, Model I has /i = 0.6, 
Qq = 0.7, Qrn = 0.3, wq = —1 and q = —0.5 today, Qq = 0.042 at nucleosynthesis. 
Similar values are found for Model II, namely h = 0.6, Qq = 0.7, Q^n = 0.3, wq = —0.9 
and q = —0.4 today, flq = 0.042 at nucleosynthesis. 

We have studied the nature of our solutions for a rather broad range of parameters of 
the potential. We have found that it is possible to obtain permanent or transient vacuum 
domination, satisfying present bounds on observable cosmological parameters, for various 
combinations of the potential parameters. 



3 Conclusions 

We have discussed a potential involving two coupled scalar fields that leads to an inter- 
esting asymptotic quintessence behaviour and accounts for data arising from high redshift 
Type la Supernovae, CMB and cluster mass distribution data. The late time dynamics 
arising from this potential, Eq (|l]), is fairly rich as it gives origin to models where ei- 
ther the acceleration regime is eternal (Model I) or instead a regime where acceleration 
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is transient (Model II) although a (non-realistic) late time non-accelerating dynamics is 
also possible. A relevant issue of this proposal is that it allows evading the conclusions of 
Refs. jT^ , |l3l , in what concerns the stability of a supersymmetric potential. Furthermore, 
since, in Model II, acceleration is transient and occurs only at present, this model is con- 
sistent with the underlying framework of string theory as is does not present cosmological 
horizons that are associated with eternally accelerating universes. 

We conclude that the late time dynamics arising from this two-field quintessence model 
is consistent with the observations as well as the theoretical requirements of stability of 
the supersymmetric ground state and the asymptotic behaviour of string theory states, 
provided the observed accelerated expansion of the universe is transient and decelerated 
expansion is soon resumed. This solution has been recently proposed to solve the con- 
tradiction between accelerated expansion and string theory [^, on general grounds; in 
this work, we have presented a concrete example of a two-field model that exhibits this 
desirable feature. 
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